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Abstract 

Recent data on B meson mixings and decays are, in general, in accord with the standard 
model expectations, except showing a few hiccups: (i) a large phase in Bg mixing, (ii) a significant 
difference (> 3.5a) between CP-asymmetries in B^ n^K^ and Bd — > n^K^ channels, and (iii) a 
larger than expected branching ratio in Bd — > tt^tt'' channel. Wc show that selective baryon number 
violating Yukawa couplings in i?-parity violating supersymmetry can reconcile all the measurements. 
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Introduction: There is still a possibility that by the time we start analyzing the LHC data, some 
indirect evidence of new physics would pop up from B meson mixings and decays. So far, most of the 
measurements in the i?-factories are in reasonably good agreement with the standard model (SM). In 
some cases, they are not, but in most such cases the uncertainties plaguing the low energy hadronic 
phenomena prevent us from making any substantial claim for new physics (NP). But, rather than 
searching for individual solutions for these discrepancies taken separately, if we seek for a collective 
solution and observe that all or most of them can be reconciled by a single NP dynamics, then that 
indeed deserves attention. Here, we focus on three such anomalies, which we call puzzles, for each of 
which a departure from the SM expectation is noticed with a reasonable statistical significance: 

(i) The Bg mixing puzzle: A model-independent test of new physics contributing to Bg mixing was 
performed with the following parametrization: 

ASMp-2if3s I 4NP„2i(</.NP_/3,) 

where Ps = aTg{-VtsV*jVcsV*^) has the value 0.018 ±0.001 in the SM. UTfit has got two solutions [l]: 
^Bsideg) = -19.9 ±5.6 , A^^ /A^^ = 0.73 ± 0.35 ; 

^Bsideg) = -68.2 ±4.9 , yl^^/^^^ = 1.87 ± 0.06 . (2) 

The SM expectation of (pBs is zero. But the above numbers show that (pSs deviates from zero by 
more than 3.7a for the first solution, while the second solution is significantly more distant from the 
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SM expectatior0- It should be noted that here the theoretical uncertainty is small, so a statistically 
significant non-zero (pB^ would constitute an unambiguous NP signal. Combining the two UTfit 
solutions, the allowed range of the mixing-induced CP-asymmetry in the Bg system is given by € 
[0.35,0.89] at 95% C.L. [2], where S^^ = sin2(|/3,| - 0bJ. 

(ii) The ttK puzzle: The observed direct CP-asymmetries in the ttK channel [3], 

acp(^d ^ vr=FK±) = -0.097 ±0.012 , acp(5^ ^ vr°A'=^) = 0.050 ± 0.025 , (3) 

imply that Aacp = acp{B^ tt^K^) - acp{Bd ir'^K^) = 0.14 ± 0.029 differs from the naive SM 
expectation of zero at 4.7cr level. In the QCD factorization approach, Aacp = 0.025 it 0.015, which 
differs from the experimental value by 3.5c7 . This is quite reliable as most of the model-dependent 
uncertainties cancel in the difference 0]. 

On the other hand, the following CP-conserving observables, as ratios of branching ratios [3] 

Rn = — =^ ^ = 1.0 ±0.07, 4) 

2 i37e[50 ^ vrOAO] + [5° ^ n^K^] 

^ ,BnB- ^^^K-HBn[B~ ^.^K~] ^^^^^^^^^ 

Bn[B+ ^ TT+AO] + [B~ ^ 7r~K^] 

are both in excellent agreement with the SM in which each of them is expected to be unity. The 
'puzzle' seems to lie in the asymmetries. 

(Hi) The vrvr puzzle: The ratio 

is in conflict with the expected relation BTl{B^ vr^vr^) >> BTl{B^ — > vr^vr^). More specifically, 
what is expected, based on different theoretical models (naive factorization [5], PQCD [6], QCDF 
[7]), is Bn{B^ TrOyrO) ~ ©(A^) Bn{B^ tt^tt^f), while what is observed is Bn{B^ ^0^°) ~ 
0(A) Bn{B^ vr^vr^). On the other hand, 

is in good agreement with the SM. 

It was shown in [5] that only a large color-suppressed tree amplitude, with other amplitudes as 
expected in the SM, can explain the rnr and vrAT data, though such a large amplitude is hard to 
extract from short-distance dynamics. We also note that large electroweak penguin (EWP) effects can 
resolve the tttt and ttK puzzles [9], but such large EWP contributions do not arise within the existing 
theoretical models. The option of suppressing the B tt'^tt" and enhancing B branching 
ratios by pumping up the charming penguins faces a serious obstacle when confronted with the ttK 
data [To]. Again, the next-to-leading order contributions in QCD factorization approaches [7] might 



^The UTfit collaboration have presented an updated estimate at ICHEP2008 (talk by M. Pierini): (/>s^ = ( — 19 ± 
7)° U (—69 ± 7)°, which shows a 2.6a discrepancy with the SM expectation. In any case, as long as this deviation from 
the SM value remains sizable, the numerical exercise leading to our conclusion holds. We thank D. Tonelli of the CDF 
Collaboration for bringing this to our notice. 
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jack up -B — > tt'^vt" branching ratio but then — > p^p^ branching ratio goes out of control. Thus, a 
collective explanation for all anomalies is hard to obtain. 

To account for the large phase in 6 — > s transition, several new physics models have already been 
proposed [Tl]. In this short paper, we show that some selective i?-parity (more specifically, baryon- 
number) violating couplings can not only provide a large phase encountered in Bg-Bg mixing but can 
also explain the vrvr and ttK riddles at the same time. 

i?-parity violating couplings: R-parity is a discrete symmetry defined as R = (— i)3B+i+2S^ where 
B, L, and S are respectively the baryon number, lepton number and spin of a particle. R equals 1 
for all SM particles and —1 for all superparticles. Unlike in the SM, conservations of B and L in 
supersymmetric models are rather ad hoc, not motivated by any deep underlying principle. However, 
such couplings are highly constrained [12]. Here, we concentrate on explicitly broken i?-violating part 
of i?-parity violation (B-RPV) only. These are contained in the superpotential, 

w = h'r,u^Dpi, (8) 

where the antisymmetry in the last two indices implies A-^-^ = —^"kj- Our selection of B-RPV couplings 
is motivated through the following chain of arguments: 

(i) First, we take only those product couplings which contribute to Bg-Bg and B^-Bd mixings via 
one-loop box diagrams. These are ^'[1^)^112 and A^23'^i2i respectively, where i corresponds to all the 
three singlet up-type flavors. 

(m) A^'^3Aj2*2, for i = 2, contributes at tree levelto b — > ccs (B^ —>■ J/'^Ks)- This is a golden channel 
for sin 2/3 measurement, yielding sin 2/3 = 0.681 ± 0.025 [3], which is slightly lower than the SM fit 
(sin2/3)fit = 0.75 it 0.04 |13pl . Now, for any i, Af23A-2\ does contaminate sin 2/? extraction any way by 
contributing to Bd-Bf^ mixing through one-loop box graphs. But, nevertheless, we refrain from using 
A213A212 to avoid any overwhelming tree level new physics imposition on the 'sin 2/3 golden channel'. 

(Hi) For a simultaneous solution of the ttK puzzle, we expect to generate a numerically meaningful 
contribution to B^ —>■ K^tt^. The corresponding quark level process b suu is triggered by ^ii'j,\i2 
for i = 1, but not for i = 3. For this reason, we consider i = 1 only as far the combination Af]^3Aj]*2 
is concerned. Regarding the other combination A^23'^i2i' again we select the i = 1 case as only this 
choice leads to 6 ^ duu {B vrvr) at the tree level. 

(iv) Thus we are left with two combinations: A'/]^3A']'J'2 and A'/23A']'2x- These consist of three inde- 
pendent couplings: X'd-^, X112 and A'/23- The strongest constraint on X'(y^ comes from n — n oscillation: 
-^113 < 0.002 — 0.1 for triq < 200 — 600 GeV [H]. On the other hand, double nucleon decay into two 
kaons puts the most stringent constraint: A'/in < 10~^^-R"2" with R = — r, the ratio between the 

hadronic and supersymmetry breaking scale. For R ~ 10"'^, the constraint is very strong: 10"''; 
while for R ~ 10~^, it gets pretty relaxed: A'/^2 ~ 1- The upper bound on A'/23 is 1.25 arising from 
the requirement of perturbative unification. 

■^Using the recent lattice measurements of the hadronic matrix elements, Bk and (^s (see Eq. (|13|l '). the authors of 
|13) have speculated a possible role of new physics to account for the difference between the fitted sin 2/3 — 0.87 ± 0.09 
(without Vub as input) and the measured value of sin2/J, which is about 2. la lower than the fitted value. 
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B-RPV contributions to observables: The product coupling A']'^3A'-|'j^2 triggers 6 — > s transition, 
while A'/23A'^21 leads to 6 ^ d transition. We define: 

h{b ^s)^ X'l.lX'l,, , h{h ^d)^ X'i^Ai ■ (9) 

These combinations contribute to Bq-Bq (q = d, s) mixing via two kinds of box diagrams, one with 
internal d'^ quark and squark and the other with u"^ quark and squark. They are given by 
{xf = m?^/m?) [15] 

= ^T^^^bMbBs, (Soi-u) + Soix,)) , (10) 



12(<?) 1927r2M? 



where 



So{x) = - -— + -——. (11) 



1 + X _|_ 2a; log x 
x)2 (1 — x) 

Above, we have assumed the relevant squarks, uji and q/j, to be mass degenerate, and we have denoted 
the common squark mass by m. 



The product coupling h{b s) also contributes at tree level to non-leptonic B decays like b 
dds and b uus, like S+ ^ if°7r+,S+ K+tt^ B^ K^7T^,Bd K+Tr-,Bs 4>Tr'^,Bs 
iT~^TT~,Bs K^K~ and their CP conjugate decayqj- Similarly, h{b — > d) provides new tree level 
contribution to different B —>■ vrvr decay modefl Thus, different decay rates receive different amount 
of SM and B-RPV contributions, and the net amplitude in each case amounts to their coherent surE0. 
The SM amplitude is calculated in the naive factorization model [5]. Considering the uncertainties 
in any such calculation, we rely on observables which are either the ratio of branching ratios or CP- 
asymmetries (in B — > ttK modes) . For the direct CP-asymmetries to proceed we need a sizable strong 
phase difference between the SM and the B-RPV amplitudes, which may be generated from final state 
interaction and rescattering. Indeed, the weak phases of the B-RPV couplings are free parameters. 
For simplicity, we have not considered the mixing between the B-RPV operators and the SM operators 
between the scale Mw and mi,. The dominant effect, which is just a multiplicative renormalization 
of the B-RPV operator, can be taken into account by interpreting the B-RPV couplings to be valid 
at the TTife scale and not at the Mw scale (thus, one should be careful in using the constraints on the 
couplings and in comparing different limits, though the numerical differences are not expected to be 
significant). 

Numerical inputs: Unless otherwise mentioned, all numbers are taken from [3]. The measured 
values of the mass differences (AMq) are 

AMd = (0.507 ± 0.005) ps"^ , AM^ = (17.77 ± O.lO(stat) ± 0.07(syst)) ps"^ . (12) 

We require sin 2/3 to he between 0.75 ± 0.04 (the SM fit value with Vub as input) and 0.681 ± 0.025 
(measured from the golden channel B^ — > J/'^Kg). 



^Contributions from lepton-number violating A'-type couplings to CP-asymmetry in B"*" -k^K channel have been 
studied in [16]. A similar study with A' couplings affecting B X^'y channel has been performed in [T^. Note that the 
B-RPV couplings we have considered in this paper would contribute to _B — > Xs7 too, but it can be kept under control. 

''interplay between Bd-Bd and Bd tt^tt" with A'-type couplings was studied in |18) . 

^It should be noted that for simplicity of our analysis we have neglected the contributions arising from ii-parity 
conserving sector in all these cases. The leading contributions from this sector to non-leptonic B decays would come at 
one-loop order, whereas the B-RPV contributions in those decays would proceed at tree level. 
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We also use the recent lattice values of the bag factors [TU] 



/i?.V^ = 281±21 MeV 



= 1.20 ±0.06 



(13) 



:S 



and the short distance factors 



0.55 



So{xt) = 2.327 ±0.044 . 



(14) 



The relevant CKM elements are |20j 



Vtd 



8.54(28) X 10 



,-3 



Vts 



40.96(61) X 10 



,-3 



7 = (75 ± 25)° 



(15) 



while the other elements are taken to be fixed at their central values. 
Results: We proceed by making two assumptions or working conditions: 

(i) The strong phase difference between the SM amplitude and the corresponding BSM amplitude 
is the same irrespective of whether it is 5 — s or 6 — > d transition. This assumption relies on flavor 
SU(3) symmetry. 

(ii) In order to calculate the amplitudes for different non-leptonic decay modes we have followed 
naive factorization approach and considered 10% uncertainty over the SM amplitudes to cover the 
different (model-dependent) non-factorizable corrections. For — > vr'^Tr^ mode we have taken this 
uncertainty to be 20%, since the SM branching ratio for this mode is Nc sensitive [5]. 

There are five parameters which we like to constrain: the magnitude of two product couplings 
(lA'i'ssA'i'sil and \X'{^lX'lu\), their weak phases ($b = Arg (A'/a^A'/ai) and $5 = Arg (X'dlX'lu)) , and 
the common strong phase difference between the NP and the SM amplitude (Ss)- We vary all of them 
simultaneously, and constrain them by requiring consistency with the observables Aacp, Rn, Rc, Rmr, 
Ra, sin 2/3, AM^, AM, and 4>b,- We also use R = Bn{B^ tt+tt")/BTZ{B° tt+K') = 0.259±0.023 
[3] to constrain those parameters. Our results are plotted in Fig. 1 and Fig. 2. Throughout our analysis 
we have taken rh = 300 GeV; a few percent variation of it will not qualitatively alter our conclusions. 
Although we varied all the parameters simultaneously, in Fig. la we projected the allowed region in 
a two-dimensional space of the magnitude (|A'/^3A'/^2l) and phase ($s) of h{b s). The red (right- 
side) patches are allowed solutions when all the five parameters pass through the filters of AM^, 
sin 2/3, Aacp, R, Rn and R^, while the blue (left-side) patches are zones allowed by AM, and (j)B^ 
only. There are small overlaps between the allowed regions from the two sets. The overlaps signify 
a common solution for all the three puzzles. With increasing statistics and with further reduction in 
theoretical uncertainties, the overlap may increase or decrease, i.e. it may or may not be possible to 
simultaneously address all the riddles with B-RPV interactions. In Fig. lb, we displayed the allowed 
zone in the plane of $5 and 6s- We note at this stage that has four sets of solutions, one in each 
quadrant, and for each such set there is an associated patch of 6s- 

Note that Rt^j^ has been deliberately kept out of the above list of constraints. If we include it, then to 
accommodate large I3TZ{B^ tt^tt^), only two sets of 6s are allowed, one in the interval (100 — > 165)° 
and the other in (195 — > 245)°. Since 6s has been assumed to be the common strong phase difference, 
its limitations of the b ^ d sector infiltrate into the 6 — > s sector as well, thus eliminating solutions 
in the second and the third quadrants. The finally allowed values of $5 lie in the range (10 — > 60)° 
and (275 340)°. Clearly, if we relax the assumption of equality of the strong phase difference (i.e. 
a common 63), solutions in all the four regions will be allowed. 
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Figure 1: (Left panel-la): The allowed zone in the plane of the magnitude of h{b — > s) and its weak phase 
(^s) "is shown. The red patches (on the right side) are scatter plots of the allowed parameters obtained by using 
AMd, sin2/3, Aacp, Ti, Rn, Rc o,nd Ra; while the blue patches (on the left side) correspond to the space allowed 
by AMs and 4>b, only. (Right panel-lb): The allowed patches in the plane of the strong phase difference (5s) 
and <i>s are displayed. 



Fig. 2a is a zoomed version of Fig. la, except that in Fig. 2a we have included all possible constraints 
at the same time. For illustration, out of the two allowed sets of <1>5, the one within the range 
(10 — > 60)° has been shown. Fig. 2b is an equivalent description replacing the magnitude and weak 
phase of h[b — > s) by those of h{b — >■ d). Note that the constraint on \h{b d)\ is one order of 
magnitude tighter than \h{b — >■ primarily because the SM prediction of the mixing is relatively 
more precise. 
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Figure 2: (Left panel-2a): Zoomed version of Fig. la, only that all constraints are now used, and focus sed in 
the first quadrant solution of ^s- (Right panel-2b): Similar to Fig. 2a, but in the space of the magnitude and 
phase of h{b d). 



Conclusions: In this paper, we wanted to solve three puzzles in B physics, namely, the large phase in 
Bg mixing, a more than 3.5 cr discrepancy between CP-asymmetries in charged and neutral B decays 
in ttK modes, and a significantly larger than expected neutral B decay in tt'^tt" channel. Here we make 
two remarks: (i) the theoretical uncertainty in the estimation of the Bg mixing phase is small and 
hence a large non-zero phase would constitute a clinching signal for new physics; (m) but, on account 
of large hadronic uncertainties associated with the ttK and tttt modes, the discrepancies observed 
in Aacp and R-n^n, though tantalizing, are not conclusive. In fact, to get rid of these theoretical 
uncertainties as much as possible, we considered the difference between CP-asymmetries and the 
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relative branching ratios. Yet, from a conservative point of view, instead of entering into a debate 
whether the discrepancies constitute 'puzzles' or 'non-puzzles', all that we wanted to emphasize in this 
paper is that if one can figure out a new dynamics beyond the SM that causes a simultaneous and 
systematic movement of all those theoretical estimates towards better consistency with experimental 
data, then that source of new physics calls for special attention. As an illustration, we advanced the 
case of explicit baryon-number violating part of supersymmetry, and we have used only two product 
couplings, constructed out of three individual ones, to explain all the data. One should keep track of 
it in the LHC data analysis, as such interactions would give lots of final state jets. 

In fact, even within the B physics context, it may be possible to infer our choices of B-RPV 
couplings (or, similar type diquark couplings) from the following observations: the coupling h{b s) 
will contaminate Bs — > K~^K~ {b — > suu at the quark level) which is used to extract 7 = Arg (V*^) 
\21\ . but it would not affect Bg DgK (b — > sen at the quark level) which is also used to determine 
7 [22]. Any statistically different measurement of 7 between these two methods will strengthen our 
hypothesis. Moreover, either of the two methods would yield 7 different from the value extracted from 
B ttK. We stress again that the falsifiability of our hypothesis, under the assumptions spelt above, 
can be judged from Fig. la by noting that the common solution zone in the parameter space arising 
from the '-B^-set' and the other data set may shrink or expand as more data accumulate. LHCb will 
definitely shed more light to these issues. 
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